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Dark matter, the CMSSM and lattice QCD 
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Recent lattice measurements have given accurate estimates of the light and strange quark conden- 
sates in the proton. We use these new results to significantly improve the dark matter predictions 
in a set of benchmark models that represent different scenarios in the constrained minimal super- 
symmetric standard model (CMSSM). Because the predicted cross sections are at least an order 
of magnitude smaller than previously suggested, our results have significant consequences for dark 
matter searches. 
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The astronomical evidence for dark matter [l], Q 
presents modern physics with some of its greatest chal- 
lenges. We need to find ways to detect it directly 
and a number of very sophisticated searches are under- 
way [1, 0, [f| , with at times confusing results @, • To 
guide those searches we need theoretical models for what 
the dark matter might be and what cross section one 
might expect it to have for scattering from hadronic mat- 
ter. The constrained minimal supersymmctric extensions 
of the Standard Model (CMSSM— see Q for a review) 
have the advantage that they are consistent with modern 
nuclear and particle physics while incorporating the com- 
pelling concepts of supersymmetry, coupling unification, 
and viable cold dark matter These models have over 
time been tuned to ensure that they are consistent with 
the latest constraints on relic abundance 2, 11 , 12, lj| , 
with the favoured dark matter candidate being a neu- 
tralino with mass of order 100 GeV and a density of a 
few per litre. 

Extensive studies of the spin-independent interaction 
of neutralinos with hadronic matter have established 
that the cross section is determined by their scalar cou- 
pling to the light and strange quark sigma commutators 
((m u + md){uu + dd)/2 and m s (ss}); see 12|, 
references therein. Given that the favoured values 



and 
of 

these quantities have been typically 50 MeV 1J] and 300 
MeV [15[ , respectively, the dark matter cross section has 
been dominated by the strange quark content of the pro- 
ton. Realizing the importance of these quantities Ellis 
et al. recently made a plea for more accurate experi- 
mental data for them [12j . While this seems unlikely in 
the short term, at least, the answer to their plea has re- 
cently been provided from an unexpected source, through 
the study of octet baryon masses as a function of quark 
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mass. Indeed, a sophisticated chiral analysis of recent 
lattice measurements has yielded surprisingly accurate 
estimates of these sigma commutators [16| . Subsequent 
lattice computations have found consistent results [17] . 

In this Letter, we use the precise new values of the 
light and strange quark sigma commutators obtained 
in Refs. [H, 17 1 to update the cross sections predicted 
within the benchmark CMSSM models that have been 
studied for sever al y ears by Ellis et al. [l2j and their 
collaborators 1(| 12 1. We show that the sizeable reduc- 
tion in the strange sigma term from the values previously 
favoured leads to a rather dramatic reduction in the ex- 
pected neutralino cross sections, with important implica- 
tions for the interpretation of experiments to search for 
dark matter. 

Sigma terms from lattice QCD. The sigma terms of 
the nucleon, for each quark flavour, are given by the 
scalar form factors evaluated at t = 0, denoted by 



m q (N\qq\N). 



(1) 



The nucleon scalar form factor is difficult to directly 
probe in experiment. The light-quark sigma term, o~e = 
(m u + rrid)(uu + dd)/2, has most reliably been accessed 
by invoking a chiral low-energy relation between tt-N 
scattering and the scalar form factor [3, EH [HI > 

(2) 



SttJV = = — Ar 



The remainder term, Ar, describes a correction to the 
low-energy theorem and is estimated to be less than 
2 MeV [3, [2(| . The shift in the scalar form factor can 
be inferred from a dispersion analysis and found to be 



rather large [21j, A CT = at{2m%) - at = 15.4 ± 0.4 MeV. 



is the Born-subtracted, isoscalar ttN scattering 
amplitude evaluated at the (unphysical) Cheng-Dashen 
point. An early experimental extraction gave = 
64 ± 8 MeV, to be compared with more a recent determi- 
nation = 79 ± 7 MeV (H. These two values lead to 



light-quark sigma terms of 



ryoungOanl . gov 



ere 



45 ± 8 MeV and 64 ± 7 MeV, 



(3) 



2 



respectively. These are to be compared with the recent 
lattice QCD determination a e = 47 ± 9MeV 0. While 
this lattice analysis tends to favour the lower value, it is 
not inconsistent with the higher extraction. We use the 
lattice determination in the current analysis. 

Extracting the strangeness sigma term is significantly 
more challenging, because the same prescription would 
lead to both a poorly converged low-energy relation and 
a large extrapolation of K-N scattering to the Cheng- 
Dashen point — see Reya A much more practical 

approach has been to resolve the pat terns of SU(3) break- 
ing among the baryon octet 1^, 25|, 2(|. In essence, the 
baryon masses give guidance with respect to the symme- 
try breaking component 



(To = mi(N\uu + dd — 2ss\N) , 



(4) 



with mg = (m u + m,d)/2. The commonly reported value 
(To = 36 ± 7 MeV is based on the phenomenological anal- 
ysis of corrections to the linear mass-splittings relations 
|25j |. with further uncertainties from yet higher order 
quantified by [26| . The difference between the estimated 
(To and the extracted ai then gives a best estimate for 
the strange-quark sigma term, as related by 



(Tf 



m s 
2ni( 



tN 



(5) 



In contrast, lattice QCD allows one to directly study 
the quark mass dependence of the baryon masses from 
first principles. New lattice results, based on an SU(3) 
baryon mass analysis [l6| and a novel application of the 
Feynman-Hellman theorem at the correlator level [l7j |. 
find cr s = 31 ± 15 MeV and 59 ± 10 MeV, respectively. 
An important feature of these calculations is the demon- 
strated internal consistency, where both studies have suc- 
cessfully shown the ability to predict the mass-splittings 
between simulations with different strange quark mass 
parameters. We argue that these results show a vast 
improvement over the earlier phenomenological analyses. 
This, in turn, implies that it is appropriate to update 
the earlier estimates of dark matter cross sections. In 
the present analysis, we assume that the reported uncer- 
tainties in the lattice QCD determinations of the strange 
quark sigma term are independent and use a naive aver- 
age of the two results, a s = 50 ± 8 MeV. We report the 
results of our analysis at the 95% confidence level to en- 
sure that the conclusions are minimally sensitive to our 
input. 

Constrained MSSM. The constrained minimal super- 
symmetric standard model (CMSSM) is inspired by su- 
pergravity mediated scenarios of spontaneous supersym- 
metry breaking, in the context of supersymmetric grand 
unified theories 0. Then one has a universal soft scalar 
mass mo, a universal gaugino mass mi/2, a universal tri- 
linear scalar coupling Aq (which is set to zero in all the 
models considered here), and a higgsino mass fi. In prac- 



tice one specifies the ratio of vacuum expectation val- 
ues tan/3 = v u /vd of the Higgs fields H u and Hd at the 
low energy scale (m^) and determines /i from this con- 
straint; the sign of \x is also a parameter, and unification 
of gauge couplings is imposed. One advantage of the 
CMSSM class of models is its simplicity, allowing for de- 
tailed scans over parameter space. The benchmark mod- 
els "A-M" of [13, El were selected to be representative 
of the regions of parameter space that yield neutralino 
dark matter with the correct abundance, consistent with 
WMAP constraints. 

Because of the incredible accuracy of the WMAP re- 
sults for relic cold dark matter abundance, a fine tuning 
of high scale parameters is typically required in order to 
have the lightest supersymmetric particle (LSP) fit this 
constraint. However, the dark matter cross section it- 
self only changes by a few percent when this tuning is 
done. For instance, the tuning that we made (described 
below) in Model L, which is a 5% change to mo, only 
results in a 2% change in the spin-independent cross sec- 
tion. Thus the direct detection predictions are insensitive 
to the changes that may have to be made to bring the 
relic abundance into line with any future updates. Thus 
the cross section results presented here are robust in that 
sense. 

Cross sections for benchmark models. We have com- 
puted the spin-independent cross section for neutralino 
dark matter for Benchmark Models A-M of [13, 11 1- 
Three of the benchmark models were also studied in 12| , 
where the reader may find the cross section formulae that 
we likewise use. While Ellis et al. use a private renormal- 
ization group evolution code, we have used the publicly 
available Softsusy [13] to compute the running parame- 
ters between the grand unified scale and the electroweak 
scale. (We did, however, compare our Softsusy results 
with those obtained using the code of Ellis et al. and 
found good agreement.) Minor modifications were neces- 
sary in two benchmark models, because of the slight dif- 
ferences between the renormalization group codes. These 
small shifts in high scale inputs were made in order to 
avoid a stau (f) LSP for models that are on the edge 
of the f exclusion region of CMSSM (i.e., dark matter 
should be neutral). In Model J we changed mo from 285 
GeV to 290 GeV. In Model L we changed m from 300 
to 315 GeV (305 and 310 were insufficient). As discussed 
above, these changes have an insignificant effect on the 
cross section values. 

Our results for the spin independent cross section with 
the proton are shown in Figs. 1-3. In all cases we see that 
there is only a mild dependence on S t tv- Furthermore, 
most models have rather small cross sections, crgi < 10~ 9 
pb. We do not show Model L because it is essentially de- 
generate with Model I. Nor do we show Model D because 
it has such small cross sections (asi < 10 -14 pb) that 
it is unobservable in the forseeable future. We also com- 
puted the cross sections for the neutron and found them 
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FIG. 1: Proton, benchmark models (labeled by letters A, 
B,...), 95% CL using lattice inputs for sigma commutators. 
Remaining models are shown in other figures below, since they 
would overlap with these. The neutron predictions are not 
shown because they are basically degenerate with the proton 
cross sections on the logarithmic scale that is shown. 



to be essentially degenerate with those of the proton. By 
way of contrast with previous estimates, in Fig. 4 we 
show <tsi for Model C using ae and a s extracted from 
lattice QCD with the determination of a s inferred from 
the phenomenological estimate for <tq = 36 MeV, using 
Eq. ((3J) above. Not only is the cross section calculated 
here, on the basis of the latest lattice QCD results, far 
more weakly dependent on H„n, but it is also typically 
much smaller. 

Details of quark flavor contributions. It is interesting 
to investigate the breakdown of how each quark flavor 
contributes to the overall cross section. First, we note 
that for the proton a si oc |/ p | 2 , where 




q—u 7 d : s 



The coefficients a^ q are determined by squark exchange 
and mixing coefficients of the LSP neutralino. The di- 
mensionless sigma commutators = o~ q /m p are taken 
from the lattice results, as explained earlier, whereas f^ G 
addresses the heavy flavor sigma commutators through 
SVZ relations, as discussed in the review [9(. We also find 
it interesting to compute the relative contribution of each 
quark flavor, defining f p = J2 q fq- F° r instance, the bot- 
tom quark contribution f£ = [2/27)m p f^ G a^ > /mb turns 
out to be rather large. Its increased importance in our 
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FIG. 2: Cross section estimates (95% CL) for benchmark 
models C, I, J. 
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FIG. 3: Cross section estimates (95% CL) for benchmark 
models A, G. 



computations can be traced to the fact that the lattice 
results give a smaller result for , the dimensionless 
strange sigma commutator. This enhances f^Q, as can 
be seen from Eq. ([6]). Results are given in Table [J for the 
point with maximim a si, which occurs in Model L. 

Conclusions. In this Letter we have shown that the 
recent lattice results [3, [13] for oi and a s have a dra- 
matic effect on predictions for direct detection of neu- 
tralino cold dark matter. Furthermore, the theoretical 
uncertainty is considerably smaller than in the traditional 
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FIG. 4: A comparison of our results (solid ellipse) for Model 
C, versus the traditional approach (dashed line) which relates 
the strange quark sigma commutator to the light quark one 
through Eq. © with <r = 36 MeV. 
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TABLE I: Example breakdown of quark flavor contributions 
to a si ■ This is for Model L at the maximum value of its cross 
section, within our 95% CL region. 



approach, as was highlighted in Fig. 4. Unfortunately, 
the most important effect of the improved estimates is the 
reduction of the strange quark sigma commutator, which 
leads to rather small dark matter cross sections. On the 
other hand, even within the 95% confidence intervals, 
our estimates come with small errors; hence any obser- 
vation of dark matter would be highly selective amongst 
the benchmark models that we have considered. In fu- 
ture work we hope to report on dark matter cross section 
predictions in models other than the CMSSM, such as 
the nonuniversal Higgs mass (NUHM) models, which are 
known [3] to allow for cross sections which are larger 
than in the CMSSM. 
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